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By extending the original Anderson singular gauge transformation for static vortices to two 
mutual flux-attaching singular gauge transformations for moving vortices, we derive an effective 
action describing the zero temperature quantum phase transition from d-wave superconductor to 
underdoped regime. In this action, quantum fluctuation generated vortices couple to quasi-particles 
by a mutual statistical interaction with statistical angle 8 = 1/2 and a dynamic Doppler shift term, 
the vortices are also interacting with each other by long-range logarithmic interactions due to charge 
fluctuation. Neglecting the charge fluctuation first, we find that the mutual statistical interaction is 
exactly marginal. In the underdoped regime, the quasi-particles are described by 2 + 1 dimensional 
QED; in the superconducting regime, they are essentially free. However, putting back the charge 
fluctuation changes the physical picture dramatically: both the dynamic Doppler shift term and 
the mutual statistical interaction become irrelevant short-ranged interactions on both sides of the 
quantum critical point. There are no spin-charge separation and no dynamic gap less gauge field in 
the Cooper-pair picture. The formalism developed at T = is applied to study thermally generated 
vortices in the vortex plasma regime near the finite temperature KT transition. The important 
effects of the AB phase scattering and the Doppler shift on ARPES data presented in Phys. Rev. 
Lett. 87, 227003 (2001) are also briefly reviewed. 



I. INTRODUCTION 

In this paper, we are trying to study the nature of 
zero temperature quantum phase transition from d-wave 
superconductor at x > x c to the underdoped regime at 
x < x c of the high temperature superconductors ( Fig.l). 





underdoped regime d-wave superconductor 

Fig 1: The temperature ( T ) versus doping ( x ) diagram of 
high T c cuprates. 

It is well known that inside the superconductor phase 
x > x c , there are low energy quasi-particle excitations 
near the four nodes of the Fermi surface, the quantum 
phase fluctuation of the Cooper pair condensates are sup- 
pressed, the positive and negative vortices are bound to- 
gether. As the doping decreases from the right to x c , the 
phase fluctuation increases. At the underdoped regime 
x < x c , the quantum phase fluctuations are so strong 
that they generate free hc/2e vortices, destroy the long 
range phase coherence of the d-wave superconductor, 
therefore the superconducting ground state. However, 
the local short-range pairing still exists, the low energy 



quasi-particles at the four nodes remain. This quantum 
phase transition is driven by the condensation of hc/2e 
vortices. Starting from this physical picture, Balents, 
Fisher and Nayak discussed Quantum phase transitions 
from d-wave superconductor to underdoped regime and 
to the anti-ferromagnet regime with focus on the con- 
densation of hc/e vortices in ffl. Earlier discussions on 
the stability of hc/e vortex starting from Mott insulating 
phase were given in Ref. [0 . The present paper focus on 
the condensation of hc/2e vortices. 

Mathematically, we can describe the transition by 
starting from the d wave BCS Hamiltonian inside the 
superconducting state x > x c in the presence of external 
magnetic field A M: 
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with di(x) = Cf(x), di(x) = c\(x) ||, V(x) is random 
chemical potential due to non-magnetic impurities which 
will be set to be zero in the following. 

Based on the earlier work on d-wave superconductor 
by Simon and Lee |3j and on t — J model by Ye and 
SachdevB, the author derived the gauge-invariant form 
in Ref. 
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where 4» is the phase of the order parameter A(r) and p x 
and p y are the momenta relative to the origin. 

Because high T c cuprates are strongly type II supercon- 
ductors ( k = A/£ » 1 ), we assume the gap amplitude to 
be a constant A(x) ~ A e^ across the transition from 
x > x c to x < x c . The phase <f> can be set to be zero 
inside the superconducting state x > x c because its fluc- 
tuation is suppressed. However in the underdoped regime 
x < x c , there are strong phase fluctuation which leads to 
the free hc/2e vortices located at Ri with parity rn = ±1, 
then the phase cf> satisfies V x V0 = 2irzJ2i niS(r — Ri). 

A lot of work has been done on a closely related prob- 
lem where quasi-particles are coupled to the vortices 
generated by external magnetic field inside the super- 
conducting state |§|,||JT^P@^,||@ . Employ- 
ing semi-classical approximation, Volovik pointed out 
that the circulating supercurrents around vortices in- 
duce Doppler energy shift to the quasi-particle spectrum, 
which leads to a finite density of states at the nodes || . 
This effect ( Volovik effect) has been employed to at- 
tempt to explain the experimental observations [fl6f of 
longitudinal thermal conductivity k xx in || . 

Starting from BCS Hamiltonian, Wang and MacDon- 
ald performed a first numerical calculation on quasi- 
particle spectrum in vortex lattice state 0. By phc- 
nomenological scaling arguments, Simon and Lee (SL) 
proposed the approximate scaling forms for longitudinal 
and transverse thermal conductivities for dirty d wave su- 
perconductors in the mixed state ||. Anderson em- 
ployed the first single-valued singular gauge transforma- 
tion to study quasi-particle dynamics in the mixed state. 
Unfortunately, Anderson made an incorrect mean field 
approximation which violates the " Time- Reversal" sym- 
metry H , therefore leads to the incorrect conclusion that 
there is Landau level quantization of energy levels of the 
quasi-particle. Franz and Tesanovic employed a different 
single- valued singular gauge transformation and studied 
the quasi-particle spectrum numerically They did 

not see the signature of Landau level quantization in their 
numerical calculations. More detailed studies in vari- 
ous kinds of vortex lattice states were conducted in Refs. 
|l2| , fl3f . The quasi-particle spectrum around a single vor- 
tex, paying special attention to the strong anisotropy of 
Fermi and gap velocities, was given in [p^ , |l5| . 

In the vortex lattice state, the vector potential only 
provides a periodic potential instead of scattering quasi- 
particles, its effect on the energy spectrum and other 
physical quantities is not transparent. In Ref. H, the 
author studied the quasi-particle transport in the dis- 
ordered vortex state where the random AB phase scat- 
tering may show its important effects. I observed that 
because infinite thin hc/2e vortices do not break the T 
reversal symmetry p), there is no Landau level quanti- 



zation, any correct mean field theory should respect this 
T symmetry (The other way around need not be neces- 
sarily true, for example, in v = 1/2 Fractional Quantum 
Hall system ^7|, the mean field state is a Fermi surface 
of composite fermion which respects T symmetry even 
the original electron system breaks the T symmetry due 
to the external magnetic field). I applied the Anderson 
singular transformation Eqn.^ to disordered vortex state 
and found that the long-range logarithmic interaction be- 
tween vortices suppresses the fluctuation of supcrfluid ve- 
locity (scalar potential), but does not affect the fluctua- 
tion of the internal gauge field. Therefore the scalar field 
acquires a " mass " determined by the vortex density, but 
the gauge field remains " massless " . The quasi-particle 
scattering from the "massless" internal gauge field domi- 
nates over those from the well-known "massive" Volovik 
effect and the non-magnetic scattering at sufficient high 
magnetic field. This dominant scattering mechanism is a 
purely quantum mechanical effects which was overlooked 
by all the previous semi-classical treatments In fact, 
it is responsible for the behaviors of both k xx and K xy in 
high magnetic field observed in the experiments [ pH| . 

When the vortices are generated by quantum fluctu- 
ations themselves are moving around, new physics may 
arise. For example, moving hc/2e vortices does break 
T reversal symmetry and leads to Hall voltage drop, a 
moving hc/2e vortex is different from a moving —hc/2e 
vortex. This in-equivalence renders FT gauge in Ref. 
]ll| useless for moving vortices although it is equiva- 
lent to Anderson's gauge for static vortices. Anderson 
gauge has to be employed to study moving vortices. In- 
tuitively, from Newton's third law, the vortices must 
also feel the counter- acting AB phase coming from the 
quasi-particles. Therefore there is a mutual AB phase 
scattering between vortices and quasi-particles. The dy- 
namic Doppler shift term on the quasi-particles can also 
be equivalently viewed as an effective gauge field cou- 
pled to the phase fluctuation (see EqnjlJ). In this pa- 
per, we enlarge Anderson singular gauge transformation 
to two mutual flux-attaching singular gauge transforma- 
tions. We then apply them to study the physics of quasi- 
particles coupled to moving vortices near the zero tem- 
perature quantum critical point. In canonical quanti- 
zation and in Coulomb gauge, we perform the two sin- 
gular transformations (which are dual to each other) 
to quasi-particles(which are spinons) and moving vor- 
tices (which are holons) respectively to satisfy all the pos- 
sible commutation relations. Just like a conventional 
singular gauge transformation leads to a conventional 
Chern-Simon (CS) term, the two mutual singular gauge 
transformations lead to a mutual CS term. This elegant 
mutual CS term describe the mutual AB phase scattering 
between vortices and quasi-particles. Alternatively, in 
path-integral presentation, the effective action describes 
the quasi-particle moving in both vector and scalar po- 
tentials due to the phase fluctuations of quantum gen- 
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erated vortices. By a duality transformation presented 
in Refs. 0,^8| to a vortex representation, the quantum 
fluctuation generated vortices couple to quasi-particles 
by a mutual CS term, the vortices are also interacting 
with each other by long range logarithmic interactions 
due to charge fluctuation described by a Maxwell term. 
Neglecting the charge fluctuation and its associated dy- 
namic Doppler shift term first, we find a fixed line charac- 
terized by the mutual statistical angle 9 and calculate the 
universal spinon, vortex and mutual Hall drag conduc- 
tivities which continuously depend on 6. This transition 
can also be viewed as a simple example of confinement 
and deconfinement transition. In the disordered phase, 
the quasi-particles are described by 2 + 1 dimensional 
QED. In the superconducting side, they are essentially 
free. We stressed explicitly the lack of periodicity of U{\) 
mutual CS gauge theory. When the U(l) charge fluctua- 
tion is taken into account, we treat both U(l) gauge fluc- 
tuations on the same footing and find the U(l) charge 
fluctuation turns the fixed line into a fixed point sep- 
arating the d-wave superconductor and some unknown 
charge ordered state. Both the dynamic Doppler shift 
term and the mutual statistical interaction become ir- 
relevant short-ranged interactions on both sides of the 
QCP. The spinon and holons are confined into electrons 
and Cooper pairs in the condensation of hc/2e vortices, 
in contrast to the condensation of double strength hc/e 
vortices discussed in JI],^| and reviewed in most general 
form in appendix B. There are no gapless gauge field 
fluctuations in the final theory, in contrast to f7(l) and 
ST/ (2) gauge theory pjj ]. We also discuss the properties 
of the two stable phases: disorder and superconducting 
phases around the critical point. 

As stressed in Ref. ||, in weakly type-II limit £ < 
A < d v , the superfluid velocity vanishes in the interior of 
the superconductor, the Doppler shift effect is completely 
absent, the AB phase scattering effect becomes the sole 
scattering mechanism Pq| . In this case, the interaction 
between vortices become short-ranged. The discussions 
in H on high T c cuprates are limited to H c \ <C H -C H c2 
where £ <C d v -C A. Except H is extremely close to 
H c i, we can safely take A — > oo limit. The interaction 
between vortices become long-ranged logarithmic inter- 
action and played leading roles as demonstrated in ||. 
In this paper, the quantum or thermal fluctuation gen- 
erated free vortex density scales as n/ ~ £~ 2 ( here £ c is 
the correlation length rather than the coherence length 
£). For high T c cuprates, k = A/£ ~ 60, £/a ~ 5, then 
A ~ 60£ ~ 300a, so as long as d v ~ £ c < A, namely ex- 
cept we get extremely close to the phase boundary, we 
can safely set A — ► oo. Then we have to consider the 
long-range interaction between the vortices which is me- 
diated by charge gauge field, its important role will be 
demonstrated explicitly in Sec. IV. All the present ex- 
perimental situations lie in the regime £ c < A. 

The paper is organized as the following. In the next 



section, we derive the effective action with both the mu- 
tual C-S interaction and the charge fluctuation, first in 
the path integral language, and then in canonical quan- 
tization representation. We stress the importance of the 
dual singular gauge transformation on the moving vor- 
tices to keep the mutual statistics invariant and derive the 
important expression of real electron operator in terms of 
spinon and holon (Eqn.^5|). In Sec. Ill, neglecting the 
charge fluctuation, we concentrate on the effect of the 
mutual C-S interaction and calculate the critical expo- 
nents and universal conductivities for general statistical 
angle 9. Implications for double layers Quantum Hall 
System are given. Some technical details are relegated 
to appendix A. In appendix B, we study double strength 
hc/e vortices in the most general form by assuming mov- 
ing hc/e vortices also scatter quasi-particles by a mutual 
C-S term with any integer statistical angle. In Sec IV, we 
discuss the combined effect of mutual C-S interaction and 
charge fluctuation, we show how the long-range charge 
fluctuation change the physical picture in Sec. Ill dra- 
matically. There are only electrons and Cooper pairs in 
the spectrum; the fixed line is replaced by a quantum crit- 
ical point (QCP) separating d-wave superconductor and 
some still unknown charge -ordered state; there are no 
spin-charge separation and no gapless gauge field fluctua- 
tions in the final theory, in contrast to U(l) or S(2) gauge 
theories investigated in Ref. pjj. We also discuss the 
properties of the two phases around the QCP. In section 
V, we apply the formalism developed for quantum gener- 
ated vortices to study thermally generated vortices and 
random static vortex array in an external magnetic field 
inside d-wave superconductor and recovered the results 
established previously in Ref. Finally, in Sec. VI, we 
discuss the connection and differences of our present ap- 
proach to the Z2 gauge theory and point out some open 
problems. 



II. THE DUAL SINGULAR GAUGE 
TRANSFORMATION, EFFECTIVE ACTION 
AND SINGLE ELECTRON OPERATOR 

In the following, we follow the notation in Ref. jl0]||] 
and introduce the spinon by performing a general single- 
valued Anderson singular gauge transformation d = Ud s : 

( e i( ^ A A 
H S = U- 1 HU, U=[ Q e _% B j (3) 

where (\>a = 0, 4>b = <f> or vice versa and <p is the phase 
of the Cooper pair. In the former (latter), the spinon is 
electron-like ( hole-like). The original Anderson singular 
gauge transformation [jl0| is devised for static vortices. 
Here we extend it to moving vortices whose phase <fi is 
also fluctuating, therefore, depends on both the space 
and time. 
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A. Path integral formulation 

Expanding H s around the node 1 where p = (pf,0), 
we obtain H s = Hi + H c . The corresponding linearized 
quasi-particle Lagrangian C qp in the presence of the ex- 
ternal gauge potential in the imaginary time t is: 

£q P = iAWr ~ ia T ) + VfiPx - a x )r 3 + v A (p y - a y )T 1 ]i/'i 
+ ^l^VfV^f) + iiplr^vrir) + {1^2,x^y) (4) 

where Vf,v A are Fermi and gap velocities respectively, 
f M = §c^0 — = x,y,r is the dynamical gauge- 

invariant superfluid momentum, it acts as a scalar scat- 
tering potential and provides dynamical Doppler shift to 
the quasi-particles (this effect can be considered as a dy- 
namic Volovik effect), = ^^((pA — <I>b) is the dynam- 
ical AB gauge field due to the phase winding of vortices. 
Note that the external gauge potential only appear ex- 
plicitly in the superfluid momentum v^. The fermion at 
node 1 (2) only couples to the a:(y)-component of super- 
fluid momentum, because the Fermi momentum at node 
1 (2) is along the x{y) direction. Because there are equal 
number of positive and negative vortices, on the average, 
the vanishing of and a M is automatically ensured in 
the Anderson gauge. 

We get the corresponding expression at node I and 2 
by changing Vf — > — Vf, v A —> —va in the above Eqn. 

£ l qp = V4[(#t ~ ia T ) - v f (p x - a x )r 3 - v A (p y - a y )r 1 ]?Ai 
- tpltf;ivfv x (f) + iiplT 3 ipiv T (r) + (1 — ► 2, x — > y) (5) 

The curvature term H c can be written as: 

H c = -[{u a ,v a } + ^^-T 3 + p-{iL c ,-n y }T 1 ] (6) 

m 2 2eF 

Where II = p + a is the covariant derivative |2j| . H c 
takes the same form for all the four nodes. Although H c 
is very important for the thermal Hall conductance in 
the presence of external magnetic field || , it is irrelevant 
near the quantum critical point, so we will not discuss it 
anymore in the present paper. 

Performing a P — H transformation ^ lct = e a pip^„ ( 
a, f3 are p-h indices ) at node 1 and the corresponding 
transformation at node 2, the above Eqn. becomes: 

£ qp = V\[{®t + ia T ) + Vf(Px + a x )r 3 + v A (p y + a y )T X ]i)i 
+ i>\$iv f v x {r) + i$T 3 ^i« T (r) + (1 -► 2, a? -» y) (7) 

In order to make the final expressions explicitly SU (2) 
invariant, we perform the singular gauge transformation 
"012a = e~ tl -^ A ~^ B ^i a and the corresponding transfor- 
mation at node 2, then a M — * — a M , Eqn.[?] takes the same 
form as Eqn.|| |23]. Adding the two equations leads to: 



£qp = i'laii^r - ia T ) + v f (p x - a x )r 3 + v A (p y - a^r 1 ]^ 
+ i>L^iaVfV x (r) + iip\ a T 3 ^ la v T {r) + (1^2,x^y) (8) 

where a = 1, 2 is the spin indices. t' s matrices are acting 
on particle-hole space. As intended, the above Eqn. is ex- 
plicitly SU(2) spin invariant. Equivalently, we can start 
with the explicitly spin SU(2) invariant approach advo- 
cated in Ref. |jj and perform the singular gauge trans- 
formation in the p-h space. 

C qp enjoys gauge symmetry U(l) u x U s (l) ( in fact, 
for static vortices, it is U(l) x ), the first being 
uniform and second being staggered gauge symmetry: 
Uniform (or external) U u {\) gauge symmetry 

<\>A -> <pA + X, <\>B -> <pB + X (9) 

Under this uniform U(l) transformation, the correspond- 
ing fields transform as: 

-» (f) + 2x, A a -> A a + d a \ 
v a — > v a , a a — > a a (10) 

d,v a ,a a all are invariant under this external U(l) 
transformation. Therefore the spinon d a is charge neutral 
to the external magnetic field. 

Staggered (or internal) U S {1) gauge symmetry 



(f>A^4>A + X, <Pb^(I>b-x (n) 

Under this internal U{1) transformation, the correspond- 
ing fields transform as: 

— > 4>; A a ^> A a 

+ d aX (12) 

Although the spinon d a is charge neutral to the ex- 
ternal magnetic field, it carries charge 1 to the internal 
gauge field a a . 

It is easy to realize that U u (l) acts only on the boson 
sector, since the spinon is charge neutral, U s (l) acts only 
on the fermion sector. In fact, U s (l) should be a dis- 
crete local Z2 symmetry for static vortices, because up 
and down static hc/2e vortices are equivalent and do not 
break T symmetry || . 

The phase fluctuation is simply 2 + 1 dimensional X-Y 
model (For simplicity, we neglect the possible anisotropy 
in the spin-stiffness): 

C ph = jvl = j(d^-2A„f (13) 

After absorbing the scalar potential scattering part 
(the dynamic Doppler shift term) into £ p h, we can write 
the total Lagrangian C = C qp + C p h as: 
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£ = ^XaliPr - ia T ) + Vf(p x - a x )r 3 + v A (p y - a a )r 1 ]V'i a 



+ {l^2,x^y) + ^{d^-A e Jff 



(14) 



where A e / f = 2A^ 



K 1 J M and the quasi-particle elec- 
tric current is : Jq = ipjr 3 ij;j, J x = Vfip\tpi, J y = vfip^i- 
It is easy to see that the dynamical Doppler shift term on 
the quasi-particles can also be equivalently thought as a 
dynamical effective gauge field coupled to the phase fluc- 
tuation. (Note that in the static vortex array generated 
by external magnetic field H, the superfluid momentum 
is static and causes Doppler shift to the quasi-particle, 
there is no feed-back effect on the superfluid momentum 
from the quasi-particles) 

From Eqn.[w|, it is easy to identify the two conserved 
Noether currents: spinon current and electric current. 
The spinon current is given by: 

Jo = ^t(x)ipi(x) + i/}l(x)ip 2 (x) 

j* = %p\(x)v F T 3 ijji(x) + V'|(a;)wAT 1 -02(a;) 

jy = tl>\{x)v A T 1 l]n{x) + tpl(x)VFT 3 l/j2(x) (15) 

Obviously the spinon current only comes from quasi- 
particle. In principle, the spinon current is not conserved 
due to scatterings between different nodes which lead 
to anomalous terms not included in Eqn.[l4|. However, 
the inter-node scatterings involve large momenta transfer 
Ki — Kj for i ^ j, so we neglect them due to momentum 
conservation in the long wave-length limit of the phase 
fluctuation |22f| . The exact conserved spin current j& is 
with a/2 inserted in the above spinon currents. 

The electric current is given by: 



31 



dC 
dA e J f 



K{d^ - A e Jf) = K{8^ - 2A^ + J„ 

(16) 



Where the first part coming from Cooper pair and the 
second from the quasi-particle. Although they are not 
separately conserved, their sum is. 

In the the electron-like Anderson gauge where et M = 
\d a <f) p4[ , the Noether current due to the the symmetry 
under <f> — > 4> + x can be written as: 



j* 

J fJ, 



dC 



1 



= K{d^-Al")-\f ll = j li 



(17) 



It is a combination of electric and spinon currents, 
therefore also conserved. 

Following Refs. |lj],[l8|], we perform a duality transfor- 
mation to EqnjlJ 



^(d^-A'/ff-i^d^ 



2 

1 ,2 

2K- 



(18) 



Where we have separated topological trivial spin-wave 
part and topological non-trivial vortex parts. 

Integrating out the spin- wave part, we get the conser- 
vation equation for the total current j* = — ij* In 
fact, as shown in the previous paragraphs, j* and j* are 
separately conserved. Therefore we can introduce spin 
and electric gauge fields by: 



3; 
3l 



ivxdva.% 
iv \d v a\ 



(19) 



We can also define the vortex current: 



fa = -^-euuxdudxcf) 
Zn 



ZuvxdvdX 



(20) 



where a" = d u (j>/2n is the vortex gauge field. 

Substituting the above expressions into Eqnfi1| we 
reach: 



(21) 



where = — \a s ^ is the total gauge field felt by the 
moving vortices. 

Using $ for the vortex operator and adding the quasi- 
particle part, we get the following effective action: 

£ = VicJSpT - iaf) + Vf(p x - o!/)t 3 + v A (p y - a^T^ipu 



+ |(0„ - ia* - za M )$| 2 + vm) + ^ap^af 



+ j/jL ~ iA e J f e^ vX d v a\ - fie ij d l a :j 



where V^(|$|) = 



1*1 



5*1*1 



(22) 



stands for the 



short range interaction between the vortices. The last 
term is due to the Berry phase in the boson represen- 
tation (l| which is a first order time derivative term. 
It can be absorbed into Aff* by redefining Afff — * 
A e Jf + ifiS/jo, it acts like an external magnetic field in 
the z ( namely fi = Q ) direction in the above vortex 
representation. There are M = 1 species of vortex and 
N = 4 species of Dirac fermion( 2 spin components at 
2 nodes at the upper half plane), the mutual statistical 
angle 9 = ±1/2. We also changed the notation by setting 
a V> _ a v ,a® = a s ,a e = a. The dynamic Doppler-shifted 
term iK^ 1 J a e^ v \d v a\ (The dynamic Doppler shift ef- 



fect) is encoded in the second to the last term in Eqn,22 
Note that this term only involves electrical gauge field. 
This effect was not written down explicitly in the Zi 
gauge theory fill . 
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Eqn|2^ again enjoys the gauge symmetry U(l) c x 
U(l) cs which is a dynamic generalization of U(l) u x U(l) s 
in Eqns.[l(],[l2|. The first one acting on the vortex sec- 
tor is the electric U(l) c ~ U(l) u gauge field, the sec- 
ond on both the vortex and fermion sectors is the U(l) cs 
mutual CS gauge field which is a dynamic generaliza- 
tion of U(l) s . The mutual Chern-Simon term enforces 
the constraints: 2n9j^ = e^xd^af^irOj^ = e^ u \d v a^. 
Physically, it means that when a quasi-particle encircles 
a vortex, it picks up a phase 2-k6. Equivalently when 
a vortex moves around a quasi-particle, it also picks 
up a phase 2ird. Although the conventional C-S term 
breaks T symmetry and has periodicity under 9^9 + 2 
the mutual C-S term does not break 



6U27 29 3CU31L 



T symmetry and does not have the periodicity under 
9 — > 9 + 1. For example, 6 = —1/2 is not equivalent 
to 9 = 1/2. 



B. Electron operator 



The general form of electron annihilation operator is: 



c a (x) = v H 



iKi-i 



Ai a -e aP e^ B e^- x ^ 20 ] (23) 



where i is the node index, f , 2 are p — h indices, a,/3 are 
spin indices (in the rest of the paper, we will stick to this 
notation ) and 4>a and 4>b are given by: 



i^dx 



>l; = - - I " <l-r 



Eqn|2^ reads: 



i=l,2 



-e a ^ 2 e- ift - s e-^ a * dx 4 2p ] (25) 

In principle, the electron Green function G(x,t) =< 
C a (x, t)C^(0,0) > can be calculated from the above 
equation. It is easy to see that the real electron in 
Eqn|2^ is invariant under the internal gauge transforma- 
tion f7(l) a . 



C. Canonical quantization approach 

In subsection A, we derived Eqn.^ by path integral 
approach, in this subsection, we rederive this equation 
by canonical quantization method and stress the impor- 
tance of the dual singular gauge transformation Eqn,29. 

Substituting Eqn.^ into Eqn.^J, we find that the sin- 
gular unitary transformation U takes the physically more 
transparent form: 









e -i<t>/2 



J J a*dx 



where, in first quantization form: 
1 



a* 



(26) 



(27) 



The first factor leads to the dynamic Doppler shift term 
which is not explicitly written in the Zi gauge theory. 
The second leads to mutual CS term with 9 = 1/2. To 
see this, in Anderson electron-like gauge, we rewrite 
in the second quantization form in terms of the vortex 
operator $ : 



d 2 rV(j){f — r)jQ(r) 



(28) 



is the 



where j'q is vortex number operator and 4>{f - 
angle the vector (r — f) makes with x axis. 

Obviously, d s still satisfies fermion commutation rela- 
tion. When the vortices are static, jg(x) = J^i 5(x — Si), 
Eqn.28 recovers the original Anderson singular gauge 
transformation and is the only necessary transformation. 
However, when vortices are generated by quantum fluc- 
tuations, they have their own dynamics. Eqn.^ alone 
can not satisfy all the possible commutation relations. 
This motivated us to introduce the following dual singu- 
lar gauge transformation to the vortex operator: 



i J a^dx 



(24) where, in second quantization form: 



a = 



(29) 



(30) 



Where j'q is given by EqnJT^. 

Obviously, $ s still satisfies boson commutation rela- 
tion. 

Most importantly, from Eqns.|26| and ^9], we can check 
that d s and <I> S commute with each other, so their mutual 
statistics is kept intact which is our original motivation. 
Note that for static vortices, the dual transformation in 
EqnJ|9| is unnecessary, because the static vortices are not 
dynamic variables. 

It is very instructive to compare the above transfor- 
mations to the well-known singular gauge transforma- 
tion leading to composite boson in v = 1/3 Fractional 
Quantum Hall state by Zhang, Hansen and Kivelson 
(ZHK) |2(J and composite fermion in v = 1/2 system 
by Halperin, Lee and Read (HLR) J^?]]: 



i I adx 

= e J 



(31) 



where, in second quantization form: 



G 



dVV^r-rXr") 



(32) 



where p(r) = ip'i/> = iplipc is the fermion number opera- 
tor. 

The crucial difference between Eqn.[n] and Eqns.26 2^ 
is that in the former, we attach electron's own 9 flux to 
itself (9 = 3 in v — 1/3 FQH state to change fermion 
to boson and 9 = 2 in v — 1/2 system to keep fermion 
statistics intact ), the constraint V x a(r) = 2ir9p(r) ( 
EqnJH ) in Coulomb gauge didi = leads to compos- 
ite fermion %p c coupled to conventional CS term which 
breaks T symmetry, has periodicity under 9 — > 9 + 2; 
however, in the latter, we attach vortex's 9 = 1/2 
flux to quasi-particle or vise versa, the two constraints 
V x a^{r) = 2ivej%(r), V x a*(r) = 27r6>jg(r) ( Eqns|||3C] 
) in Coulomb gauges diaf = didf = naturally leads to 
quasi-particles and vortices are coupled by mutual CS 
term in additional to the electrical gauge field which me- 
diates the logarithmic interaction between vortices. The 
mutual CS term does not break T symmetry and does 
not have the periodicity under 6—^6 + 1. 

It is easy to see that when a vortex is moving around 
a closed loop, it pick up two phases, one is U(l) cs phase 
due to the spinon current, described by the mutual C- 
S term, another is U(l) c phase due to the total electric 
charge current described by the Maxwell term. In the 
following section, we will neglect the electric charge fluc- 
tuation (therefore its associated Dynamic Doppler shift 
term) and concentrate on the mutual statistics term. In 
Sec. IV, we will study the combined effect of charge fluc- 
tuation described by the Maxwell term and the mutual 
C-S term. 

For simplicity, we neglect all the possible anisotropy in 
the quasi-particle and vortex velocities. We also take the 
relativistic form for both fermion and boson, because the 
difference between the velocity of spinon and that of vor- 
tex in Eqn.^2] is expected to be irrelevant near the zero 
temperature QCP M. 



III. THE EFFECT OF MUTUAL STATISTICS 

In this section, we neglect the charge fluctuation, 
namely setting a M = in Eqn.|2^. The associated dy- 
namic Doppler shift term also drops out. As pointed out 
in H and re-emphasized in the introduction, this is the 
weakly type II case. So far, weakly type II supercon- 
ductors are all s-wave superconductors where there is no 
low energy quasi-particles. So the results in this section 
can only be applied to presumed " weakly type II d-wave 
superconductor " which has not been discovered yet in 
nature. The model is also a relativistic version of the 
Hall drag problem in double layer Quantum Hall system 
p6[ , therefore the model itself is interesting on its own 
right and deserves detailed investigation. 



The charge fluctuation can be suppressed by assuming 
the condensation of hc/e vortices reviewed in appendix 
B. We could add the kinetic and potential terms for the 
hc/e vortex operator $2 to Eqn.E2t 



£* 2 = \(d^-i2a^ 2 \ 2 + V(\^ 2 \) 



(33) 



As shown in (lj and in appendix B, there is no mutual sta- 
tistical interaction between spinon and hc/e vortices $2- 
The long range logarithmic interaction between hc/2e 
vortices $ and hc/e vortices $2 is mediated by the elec- 
trical gauge field a M . Condensing < $2 >= $20 will 
generate a mass term % a (a^)t which dominates over the 
Maxwell term. Integrating out leads to 



J_rl 



{f£") 2 -iA e ' f e M vxB v ji + (fi)*] 



(34) 



All the generated terms only renormalize the short 
range interactions already included in V(|$|). 



A. Quantum Critical point 

In order to calculate conductivities, we add two source 
fields and A* for the quasi-particles and vortex re- 
spectively: 



C = ^H(d» - ia% - iA^ a 



\(du 



i A*)$| 2 + y(|$|) + — afe^af (35) 



where a = 1,2,3,4 stands for N = 4 species of Dirac 
fermion. In fact, two Maxwell terms for eft and a* can 
be added to the above equation, but they are subleading 
to the mutual CS term in the low energy limit. It is ex- 
pected that there is no periodicity under 9 — > 9 + 1 in 
the continuum limit. 

The RG calculation in Refs. |2^,^,|l| can be used to 
show that 9 is exactly marginal, therefore there is a line 
of fixed points determined by the mutual statistical angle 
9. In order to calculate the spin conductivity along this 
fixed line, a source field could be introduced to couple to 
the spin current = ^al^i^a^fS- Similar calculations 
follow. 

Integrating out both fermion and boson leads to: 

C = -\ai{-k)Ut„(k)ai{k) - \at{-k)Ult{k)at{k) 

ia£(-*)Il£,(fc)a?(fc) ~ \al(-k)Ult(k)at{k) 

- AU-k))e,Mat(k) - At(k)) 

^«(~ k ) A^(-k))e^ x k x (at(k) - Af(k)) (36) 

where the exact forms of n' s are dictated by gauge in- 
variance and Furry's theorem : 
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kirhj. 



A: 2 



electric gauge field a M , but noi the statistical gauge field 



n^(fc) = n^(fc) = n 2 (fc) e ^ A fc, 
n* !/ (fc) = n 3 (fc)fc(v-^) 



(37) 



Where ITi , II2 , II3 are the polarizations for 



a* a*, a* a*, a* a*. 



If we are only interested in DC conductivities, for sim- 
plicity, we can put fc = 0, Eqn.^6| becomes ( for the most 
general form, see appendix A): 



£ = --af(-w n )Ui\w n \af(<jJ n ) - -a^ (-u n )Tl2eijU> n a® {uj n ) 
I A J 

1 rfc / \ -■ — r 1 1 ffw \ 1 



- -af (-cj n )TL 3 \w n \af - -af (-w n )n 2 eyW n anu„) 

z Z J 

- 2g (af ~ Af (-w„))e i3 cj„(a*(w„) - A*(cj„)) 

- ^(o?("«n) ~ (-^))e,^„(af ( W „) - (38) 
If we define af (w„) = e y -a*(w„), Af (w„) = e y A*(w„), 



The equations [4l| , [42| are exact, but III, II2, II3 can only 
be calculated perturbatively in the coupling constant 
g 2 = 2tt9. The renormalized propagators for and a* 
can be found from the following Feymann diagrams: 

r -v O 



Fig 2: The renormalized propagators for eft (thick wiggle line ) 
and a* (thick dashed line). The thin wiggle line stands for eft , the 
dashed thin line for a* , the thin solid line for fermion propagators, 
the thick solid line for the boson propagators. 



By using the bare propagators < 
-£fiv\k\/k 2 and the bare fermion and boson loop results 



afaf, >= 



n/o 



■j^k(S Pil , — k^k u /k 2 ), we can find easily 



the above equation becomes diagonal in the spatial in- 
dices i = 1,2. Finally, integrating out af, af leads to 



the renormalized a^ and a* propagators G^ v = G* y = 
— lisjK'W — k^k u /k 2 ). In contrast to the conventional 



1 



(Af,if) 



Af 
Af 



(39) 



Where spinon, vortex and mutual Hall drag conduc- 
tivities are: 



CS theory studied in 29|3(|, the propagators are even in 
fc, this is because the theory respects T symmetry. On 
the other hand, in contrast to the Maxwell propagators, 
they behave as 1/fc instead of 1/fc 2 . 

The three loop diagrams for II 1 are given by: 



a* 


= (-? 


a* 


=<?> 2 


<j h 





III 


ITII3 


+ (1/0- 


n 2 ) 2 




n 3 




ITII3 


+ (1/0- 


n 2 ) 2 


n 2 - 


oin^-s - 


nn|) 


111II3 


+ (i/0- 


n 2 ) 2 



VVW ! 





(40) 



In fact, all the three conductivities can be written in 
the elegant connection formula: 



Pij - (PFB)ij 



(41) 



with the conductivity tensor of fermion and boson given 
by 



Fig 3: The three loop diagrams of Eli. The thick wiggle line 
stands for the renormalized propagators of aft, the thin solid line 
stands for the fermion propagator. The one loop diagram is not 
shown. 

By using G^ v and extracting the symmetric part of the 
gauge propagator in the large N results in Refs. p9| , ^0[ , 
we are able to calculate the above three loop diagrams. 
Furry's theorem can be used to eliminate large number 
of null diagrams. We get the following series: 



ofb = 



it -n 2 
n 2 n 3 



(42) 



(43) 



Although its form is similar to the conventional con- 
nection formulas discussed in |27 2j],|3(J, the physical in- 
terpretations of the conductivities are quite different (see 
the following). 

When the vortices are generated by external magnetic 
field and pinned by impurities as discussed in the 
total conductivity is the same as the fermion conduc- 
tivity because the static vortices do not contribute. As 
explained in section V, the static vortices only feel the 



where N — 4 ( M — 1 ) is due to the sum over 4 (1) 
species of Dirac fermions ( bosons ) . 

The three loop diagrams for II3 are given by: 



o 
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Fig 4: The three loop diagrams of II3. The thick dashed line 
stands for the renormalized propagator of a*, the thick solid line 
stands for the boson propagator. The one loop diagram is not 
shown. 

For bosons, only one loop result is known: 



n. 



MN-(l + g 4 + g a 



(44) 



Fig 6: The two loop diagrams for the self energies of spinon (a) 
and vortex (b). The renormalized operators are used. 

The results are: 



Note that both III and II3 are even functions of 9. 

From Furry's theorem, one of the first non-vanishing 
diagram for II2 is: 





48tt 2 



-M 



12tt 2 



(45) 



where tjxy ~ 0.038 is the anomalous dimension for the 
3d XY model ||. 

The correlation length exponent can also be calculated 
to two loops by the insertion of the operator <&^<i>: 



Fig 5: One of the four loop diagrams of 112. the thin solid line 
stands for the fermion propagator, the thick solid line for the boson 
propagator. All the other four loop diagrams can be obtained by 
shuffling the positions of the three bare propagator lines. 



MN(g 6 
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which is an odd 



The series is IT2 
function of 9. 

From EqnJfO], <7^,£r* are even functions of 6, but the 
mutual Hall drag conductivity is an odd function of 9. 
These are expected from P — H transformation. Under 
the P — H transformation of the vortex operator $ — > $t 
in Eqn.|35|, it can be shown that A* , 9 is equivalent to 
-A*, -9. From Eqnj3g we reach the same conclusions. 
Specifically, cr H takes opposite values for 9 = ±1/2, the 
periodicity under 9 — > 9 + 1 is not preserved. Exper- 
imentally, the Hall drag conductivity can be detected 
by measuring the transverse voltage drop ( or transverse 
temperature drop for thermal conductivity ) of spinons 
due to the longitudinal driving of vortices. The Hall 
drag conductivity in double layer Quantum Hall systems 
has been investigated by several authors ||(|. In double 
layer systems, the electrons in different layers are treated 
as two different species. There is a mutual CS interac- 
tion between the two species (both are fermions) which 
is directly responsible for this Hall drag conductivity, al- 
though the Coulomb interaction between the two species 
is responsible for the Coulomb ( longitudinal ) drag |J7| . 
This example shows that no external magnetic field is 
needed to produce a Hall effect ! For example 9 = ±1/2 
lead to opposite Hall drag conductivities although static 
a = ±1/2 vortex leads to no Hall conductivity. 

Again, by using the large N result of |2^,|3(|, we find 
the anomalous dimensions of the fermion and vortex to 
two loops: 



Fig 7: The cross stands for the operator insertion of <3> t <j> 

The result is: 

,4 



V = V X Y 



9 



12tt 2 



tN 



(46) 



where vxy ~ 0.672 is correlation length exponent for the 
3d XY model |33l. 



It is instructive to go to dual representation of Eqn, 35 
namely go to the boson representation: 

+ \{ft»? + ^%^(at - af/9) ~ iA% vX d v a1 (47) 

where V{\cf>\) = ra\\4>\ 2 + <? |0| 4 + • • ■, 

Integrating out a* leads to the constraint up to a pure 
gauge: 



at = 9at 



(48) 



Substituting the above constraint to Eqn.[l7] and set- 



ting = a^, we find: 



C = ^(d,-t9a,- l A^ a + -f^ 

+ - ta^l 2 + V(\ct>\) - iA*€^, x d u ax (49) 

The above Eqn. indicates that fermions and bosons 
are coupled to the same gauge field whose dynamics is 
described by Maxwell term instead of C-S term. This 
Eqn. is simply 2 + 1 dimensional combination of spinor 
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QED and scalar QED. This action is similar to a rela- 
tivistic analogue of the U(l) gauge field theory investi- 
gated extensively in Ref. EJ . The RG analysis at 4 — e 
by dimensional regularization is possible, because the 
marginal dimension of all the relevant couplings are 4. 
However a RG analysis directly at 2 + 1 dimension is 
formidable. The physical meaning of 6 is obscure and 
the exact marginality of 9 is a highly non-trival result in 
the boson representation, but all these become evident in 
the dual vortex representation Eqn.^. It is evident that 
there is no periodicity under 9 — > 9 + 1 in Eqnj49|. 

If we perform duality transformation again on Eqnj49| 
to go to the vortex representation, then we recover 
EqnJH upon neglecting the two Maxwell terms which 
are subleading to the mutual C-S term. 

In the next subsection, setting the two source terms 
vanishing, we look at the properties of the different 
phases on the two sides of this quantum critical point. 



B. 2 + 1 dimensional QED in the underdoped phase 
and free fermions in the superconducting phase 



(52) 



where m$ is the mass of the vortex. 

The Dirac fermions become free. In fact, we reach the 
same description from the boson representation Eqn, 49. 
Because in superconductor phase, the boson <f> condense, 

,2 

therefore generates a mass term -^(a^) 2 which renders 
the Maxwell term ineffective. We reach the same con- 
clusion from both sides by identifying m$ ~ 0q which is 
dual to the relation in the disordered phase ~ $q. 

In short, in the disordered phase, the system is de- 
scribed by spinor QED Eqn.|5l|; in the superconductor 
side, by free Dirac fermion. We can view the transi- 
tion as a simple example of confinement-deconfinement 
transition. In the confined ( disordered ) phase, the bo- 
son and fermion are confined together by the fluctuating 
gauge field. In the deconfined ( superconductor) phase, 
the boson condensed, the gauge field becomes massive, 
the fermion becomes free. There is a line of fixed point 
governed by the mutual statistical angle 6 separating the 
two phases. 



In the disordered phase, the vortex condense < $ >= 
$o which generates a mass term for a* in Eqn.|35| 

f («£)? (50) 

where the subscript t means transverse projection. 

Integrating out the massive a* leads to a Maxwell term 
for afr. 

C = iil,{d,~ia^ a + ^ 1 U^f (51) 

This is simply 2 + 1 dimensional spinor QED which was 
studied by large N expansion by the authors in . They 
used four component Dirac fermions and discussed possi- 
ble dynamic mass generation of the fermions which break 
the chiral symmetry, but preserves Parity and Time re- 
versal symmetry. However, for 2+1 dimensional compact 
QED, the possible important instanton effects were not 
well understood yet. The possible connection between 
this dynamic mass generation and Anti-ferromagnetism 
in the context of gapless flux phase was discussed by 
Marston Q. 

In fact, we reach the same description from the boson 
representation Eqn.[l9| Because in the disordered phase, 
the boson <f> is massive, therefore can be integrated out, 
it generates the Maxwell term — !— f?„. which dominates 
over the existing non-critical Maxwell term. We reach 
Eqnj^ after identifying ~ $q. 

In the superconductor phase, the vortex <!> is mas- 
sive, therefore can be integrated out, it leads to the old 
Maxwell term j^(f^) 2 - Integrating out a* generates 
a mass term: 



IV. THE EFFECT OF CHARGE FLUCTUATION 
AND THE ABSENCE OF MASSLESS DYNAMIC 
GAUGE FIELD AT T = 

In this section, we try to investigate the effect of charge 
fluctuation on the fixed line characterized by the statis- 
tical angle 9 discovered in the last section by considering 
the combined effects of U(l) cs mutual statistical gauge 
fluctuation and U(l) c electrical gauge fluctuation and 
treat both U(l) gauge fields on the equal footing. We 
find the condensation of hc/2e vortex condensation in- 
deed leads to the confinement of spinon and chargon into 
Cooper pairs and electrons, in contrast to the conden- 
sation of hc/e vortex. The fixed line in the last section 
is destroyed and replaced by a quantum critical point 
separating the superconducting state and some sort of 
charge ordered state. There is no gapless gauge field left 
in the final theory, in sharp contrast to the U(l) or SU (2) 
gauge theory Both mutual C-S interaction term and 
the dynamic Doppler-shift term become irrelevant short- 
ranged interactions. 

The authors in Ref. [^0| studied the similar model. 
They performed a perturbative RG calculation and 
reached very different conclusions from ours. 

A. Quantum critical point 

Putting a* — > a* — in Eqn.|2^, we get 

£ = ^(fy - iajftya + \(9, O^l 2 + V(\<f\) 

+ 2^a*e ptfA ft,aJ - ia^xMaf/e + A e x ff ) + I/*, (53) 
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Integrating out the electric gauge field a M leads to: 

£ = ^7m@m " itflWa + ~ ^*)*| 2 + 7(|*|) 
+ ^<e M , A 9,a! + i(ojf/fl + A?*)* (54) 

Comparing to Eqn.pq, it is easy to see that the charge 
fluctuation leads to a mass term for the gauge field at- 
Shifting + AffJ — * a M and adding the gauge fixing 
term ^(c^a^) 2 , we can integrate out the massive gauge 
field a M in Lorenz gauge a — and find: 

£ = iphAd, - iOA'ff)i> a + \(d, - za*)$| 2 + y(|<&|) 



(55) 



Note the Maxwell term for a* is generated by the in- 
tegration over the massive a^. 

Setting = 0, integrating out the fermions only leads 
to higher derivative terms than the Maxwell term: 



C = |(^-iaj)*| 2 + + \(fp 2 i CO! 



where • • • means higher than second order derivatives. 
Therefore the vortex and fermion are asymptotically de- 
coupled. It indicates that the charge fluctuation ne- 
glected in the last section destroy the fixed line char- 
acterized by 9. However, very different conclusions are 
reached in Ref. |BQ| . We think that the perturbative RG 
calculation in pOfl may not treat the charge gauge field 
fluctuation correctly. 

Just like the last section, it is instructive to go to dual 
representation of Eqn.|54|, namely go to the boson repre- 
sentation: 

£ = ^(fy ~ iatWa + |(0 P - iafW + 7(|0|) 



+ l Ut v ? + Ke tlvX d v {a{-ai/e) 



(57) 



Integrating out a* leads to the same constraint as 
Eqn|48 up to a pure gauge: 



4 = eat 



(58) 



Substituting the above constraint to Eqn|57| and set 



ting = dp, we find: 



£ = ^7^ - iOa^a + |(0 M - ia^\ 2 + 7(|0|) 



(59) 



Comparing to Eqn,49, the important difference is that 
the gauge field acquires a mass due to the charge fluctu- 
ation which renders the Maxwell term ineffective. Note 
that the dynamic Doppler shift term is encoded in A e J? . 



Up to irrelevant couplings, we can safely set a M = A e Jf 
in the above equation and find 

£ = i>ll»{d, ~ i26A^ a + |(a„ - i2A^\ 2 + 7(|0|) + • • 

(60) 

where • • • means the irrelevant couplings between bosons 
and fermions. 

The above Eqn. leads to the conclusion that the 
mutual statistical interaction between hc/2e vortex and 
spinons leads to the confinement of spinon and char- 
gon into electron and Cooper pair. The quasi-particles 
carry charge 29e which can be identified as electrons or 
holes, the Cooper-pairs are described by 3D X-Y model. 
The mutual C-S coupling and the dynamic Doppler-shift 
coupling between Cooper pairs and electrons are short- 
ranged interactions which are irrelevant near the quan- 
tum critical point. There is no gapless gauge fields in 
the final effective action ! However, there are several im- 
portant factors left out in the above analysis: (1) The 
Berry phase term for the boson which is first order time 
derivative term is not included. This term is vanishing 
only at half-filling with particle-hole symmetry hopping. 
(2) The long-range Coulomb interaction between Cooper 
pairs and electrons is not included. When the both ef- 
fects are taken into account, the final ground state in the 
underdoped side at T — maybe an insulating state with 
some kind of charge order. The nature of the QCP be- 
tween the d-wave superconducting state and this charge 
ordered state is an important unsolved problem. How- 
ever, the conclusion that there is no spin-charge separa- 
tion and no gapless dynamic gauge fluctuation remains 
robust ! 



B. Disordered and superconducting phases 

We follow the discussions in the previous section. In 
the disordered phase, the vortex condense < * >= $0 
which generates a mass term for a* in Eqn.[54|. Integrat- 
ing out the massive a* leads to a Maxwell term for a^. 



1 

4*| 



(61) 



Because of the mass term, the Dirac fermions become free 
and carry spin 1/2 and charge We. In fact, we reach the 
same conclusion from the boson representation Eqn.|59|. 

In the superconductor phase, the vortex * is mas- 
sive, therefore can be integrated out, it leads to the old 
Maxwell term 4^-(/*„) 2 - Integrating out a* generates 
a mass term for af,: 

£ = ^7m(^ - *4Wa + ^(af) 2 + \{a+/e + Af/f) 2 

(62) 
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where m$ is the mass of the vortex. 

Diagonizing the last two mass terms leads to a con- 
tinuously changing charge . This is expected, because 
in the superconducting state, the quasi-particle is a lin- 
ear combination of electrons and holes, therefore carry- 
ing continuously changing charge. In fact, we reach the 
same description from the boson representation Eqn. 59 
Because in superconductor phase, the boson condense 
< (j) >= 4>o i therefore generates a mass term for a M : 



C 



(63) 



Which is essentially the same as Eqnj62| after the identi- 
fication m$ ~ </>q. 

In short, in the disordered phase, the quasi-particles 
carry charge 26e and spin 1/2 which can be identified as 
electrons or holes. When the Berry phase term and the 
long-range Coulomb interactions are taken into account, 
the electrons may form some charge ordered state. In the 
superconducting side, the quasi-particles carry continu- 
ously changing charge which is a linear combination of 
electrons and holes. 



V. APPLICATION TO STATIC DISORDERED 
VORTEX ARRAY AND THERMALLY 
GENERATED VORTICES 

So far, we discussed the vortices generated by quantum 
fluctuations in the underdoped regime. Vortices can also 
be generated by an external magnetic field inside super- 
conductors or by thermal fluctuations in the pse udo-gap 
regime. They were discussed in Refs. || and pp|] respec- 
tively. In this section, we apply the formalism developed 
for quantum generated vortices at T = to study vor- 
tices generated by an external magnetic field or by ther- 
mal fluctuations. By studying the three different kinds 
of vortices in a unified picture, we not only recover the 
previous established results, but also provide additional 
physical insights on the internal gauge field fluctuations. 

Inside the d-wave superconductor, although the 
Cooper-pair condensates dominates the electrical trans- 
port at zero frequency, the quasi-particles near the four 
nodes of a d-wave superconductor are responsible for the 
heat and spin transport at low temperature. Recently 
Yc studied the quasi-particle transport in a random vor- 
tex array with a fully quantum mechanical approach Q . 
Although the previous semi-classical approach only cap- 
ture the physics of the Volovik effect, fully quantum me- 
chanical approach treats the Volovik effect ( scalar field 
scattering ) and the AB phase scattering ( the gauge 
field scattering ) on the equal footing. Ye found that 
the long range logarithmic interaction between vortices 
suppress the superfLuid velocity fluctuation, but does not 



affect the internal gauge field fluctuation. He concluded 
that the quasi-particle scattering from the random gauge 
field dominate over that from the superfluid velocity ( 
the Volovik effect). 

Roughly speaking, the mixed state with vortex density 
n v corresponds to the vortex condensed phase (disordered 
phase). As discussed in the previous sections, the vor- 
tex condensation leads to a mass term ^«(a*) 2 for the 
statistical gauge field a* in Eqn]2|. Adding the mass 
term to the equation, integrating out a* and keeping 
only the time component of the charge gauge field ao 
which mediates the long-range density-density logarith- 
mic interaction between the static vortices, we get two 
static Maxwell terms for a% , v a respectively and a mass 
term for v a : 



P 



■(/^) 3 + f(«a) 2 



(64) 



V_ 

4n v wap/ ' An v 
where p is added to stand for an unknown function of pin- 
ning potential, = d a a^ - dpa%, f^p = d a vp - dpv a . 

In Eqnj64[ the averages over a% , v a should be under- 
stood as quenched instead of annealed averages. Namely 
there is no feedback effects on the gauge field and scalar 
field propagators from the fermions f4l| ] . The above equa- 
tion shows that the long-range logarithmic interaction 
between vortices suppress the fluctuation of the super- 
fluid velocity and leads to a mass term for the superfluid 
velocity and the incompressibility of the vortex system. 
The quasi-particle is moving in a long-range correlated 
random magnetic field and short-range correlated scaler 
potential. The distribution of the random gauge field is 
given by the corresponding static Maxwell term and that 
of the scalar field is given by the corresponding static 
Maxwell term and the mass term. Eqn. |6J is essentially 
the same as Eqn. 9 in Ref. p|. 

Strictly speaking, the gauge field from the random vor- 
tex array is a random Zi gauge field. In Ref. ||, in 
order to perform an analytic calculation in the contin- 
uum limit, Ye made an approximation by replacing the 
discrete random AB phase by a continuous random 
C/(l) AB phase and performed a standard diagrammatic 
perturbation theory. 

In Ref. [Q, assuming the finite temperature phase 
transition from d-wave superconductor to pseudo-gap 
regime is a Kosterlize and Thouless (KT) vortex and anti- 
vortex unbinding transition |m| , Ye studied the electron 
spectral function in the vortex plasma regime in Fig.l 
[ fjO[ . Starting from Eqn.p5|, he found that the random 
gauge field missed in the previous semiclassical ap- 
proaches [ f42| destroyed the coherent spinon motion and 
leads to branch cut singularities and pertinent non-Fermi 
liquid behaviors. There are three subtle differences be- 
tween thermally generated vortices from the static vor- 
tices generated by an external magnetic field inside the 
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d-wave superconductor: (1) The vortex density n v should 
be replaced by free vortex density n f above KT transition 
temperature. (2) Thermally generated vortices are mov- 
ing around, the U(l) character is stronger than Z 2 char- 
acter. They destroy the superconductivity. But static 
vortices only destroy the superconductivity inside the 
vortex cores. (3) The distribution of thermally gener- 
ated vortices is precisely given by vortex hydrodynamics 
in the long- wavelength limit. The free vortex density n/ 
is the only relevant parameter. We can set p equal to 1 
in Eqng 

After taking into account these subtle differences, the 
above procedures leading to Eqn.^34| can be equally ap- 
plied to thermally generated free vortices with density n / . 
The following alternative derivation of Eqnj6^ may be 
stimulating: the vortices being treated classically, their 
commutation relations can be neglected, the dual singu- 
lar gauge transformation in Eqn.^9] on the vortices is not 
necessary, the Berry phase term for the boson ( the linear 
time derivative term) can also be neglected. More specif- 
ically, only the time component of the charge gauge field 
ao in Eqn.|2^ is kept to mediate the long-range density- 
density logarithmic interaction between the vortices. Ob- 
viously, ao couples to the vortices the same way as the su- 
perfluid velocity v a couples to the spinon, its fluctuation 
leads to a mass term for the v — v correlation. However, 
being orthogonal to the spatial component eft , oft — oft 
correlator remains gapless. This also leads to exactly the 



same equation as Eqn.64 with p setting equal to 1 and 



the same conclusions as those first reached in Ref. M . 

However, at zero temperature, the vortices being 
treated quantum mechanically, all the components of the 
charge gauge field should be kept to mediate the long- 
range current- current logarithmic interactions between 
the vortices. Then the charge gauge field fluctuation 
leads to mass terms not only for w, but also for the spa- 
tial component eft . There is no gapless dynamic gauge 
fluctuations in the Cooper-pair picture as demonstrated 
in Eqnj6^. The phenomenon that the two gapless gauge 
fields cancel each other is similar to the standard Higgs 
mechanism. In Higgs mechanism, the gapless gauge field 
is eaten by gapless Goldstone modes. The two gapless 
modes annihilate each other, the Goldstone mode disap- 
pears, while the gapless gauge field becomes massive. 



VI. DISCUSSIONS AND CONCLUSIONS 

Based on the earlier work by Balents, Fisher and Nayak 
, Senthil and Fisher developed Z 2 gauge theory [ pj|j20| 
to study quasi-particles coupled to vortices generated by 
quantum fluctuations. By breaking electrons and Cooper 
pairs into smaller constitutes: chargons with charge e, 
spin and spinons with charge 0,spin 1/2, SF intro- 
duced a local Z 2 gauge degree of freedom to constrain 
the Hilbert space to be the original one. The effective 



action describes both chargons and spinons coupled to 
local fluctuating Z 2 gauge theory with a doping depen- 
dent Berry phase term. By the combination of standard 
duality transformation of 3 dimensional XY model and 
that of Z 2 gauge theory, the action is mapped into a 
dual vortex representation where the hc/2e vortices and 
spinons are coupled by a mutual Z 2 CS gauge theory. 
As usual vortices in XY model, the hc/2e vortices also 
couple to a fluctuating U(l) gauge field which mediates 
the long-range logarithmic interaction between the vor- 
tices. Starting from the dual representation, the authors 
in Ref. |^0| studied a transition from d-wave supercon- 
ductor to confined Mott insulator driven by the condens- 
ing of hc/2e vortices at half filling. In order to study the 
critical behaviors of this particular confinement and de- 
confinement transition, they replaced the Z 2 mutual CS 
theory on the lattice by U(l) mutual CS theory in the 
continuum and performed renormalization Group (RG) 
analysis. Some of their RG analysis may not be correct 
as demonstrated in Sec. IV. In contrast to U(l) gauge 
theory which has only confined phase, Z 2 gauge theory 
has both confined and deconfined phases. The topolog- 
ical excitation of the Z 2 gauge field play an important 
role. If hc/2e vortices condense, the visons are also con- 
densed, the system is in the confined phase, the char- 
gon and spinon are confined into electron, Cooper-pair 
and magnon. If hc/e vortices condense, the visons are 
gapped, the system is in the deconfined phase, the spinon 
and chargon are asymptotically decoupled in the long- 
wavelength limit. SF argued that the cuprates are in the 
deconfined phase, therefore the zero temperature tran- 
sition at x = x c is driven by the condensation of hc/e 
vortex. They further proposed the " vison trapping ex- 
periment " to test their spin-charge separation scenario. 
Several groups performed such experiments, but found 
no signature of visons . 

The original Anderson singular gauge transformation 
ifiof was proposed for static vortices. In this paper, we ex- 
tend the one singular gauge transformation for static vor- 
tices to two mutual flux-attaching singular gauge trans- 
formations for moving vortices generated by quantum 
fluctuations. By making a close analogy to the con- 
ventional singular gauge transformation of FQH system, 
we perform the two singular gauge transformations at- 
taching flux of moving vortices to quasi-particles or vice 
versa. Just like conventional singular gauge transfor- 
mation leads to conventional CS term, the two mutual 
singular gauge transformations lead to mutual CS term. 
In this way, we propose an intuitive and physical trans- 
parent approach to bring out explicitly the underlying 
physics associated with the condensation of the hc/2e 
vortices. When considering both mutual C-S interaction 
and long-range logarithmic interaction between the vor- 
tices, we find that there are only electrons and Cooper 
pairs in the spectrum and there are no gapless gauge 
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field fluctuation in contrast to the U(l) or SU{2) gauge 
theory plj . When taking into account both the linear 
time derivative Berry phase term for the bosons and the 
long-range Coulomb interactions between Cooper pairs 
and electrons, the true ground state in the underdoped 
side may be some unknown charge ordered state. Under- 
standing the zero temperature transition from d-wave su- 
perconductor to this unknown charge ordered state is an 
important future research direction. The electron spec- 
tral weight starting from Eqn.|2^ measured by ARPES 
in the pseudo-gap regime at T > T c indicated in Fig.l 
was presented in a separate publication [Q. The the- 
ory presented in this paper focus on the confined phase. 
The recent experiment []38f indicates that the cuprates 
may be inside the confined phase. In principle, our the- 
ory should be equivalent to the Z2 gauge theory in the 
confined phase side after tracing out the fluctuating Z2 
gauge field except the Doppler shift effect missing in Z2 
gauge field is also explicitly incorporated into the present 
approach. 

The importance of the Volovik effect and the AB phase 
scattering in random static vortex array was stressed in 
Ref. 0. As concluded in this paper, when the vortices 
are generated by quantum fluctuations at T = 0, the dy- 
namic Doppler shift effect and the mutual C-S phase scat- 
tering both become irrelevant short-ranged interactions. 
However, there is a crossover regime between these two 
limits which is the vortex plasma regime near the finite 
temperature KT transition. In this regime, there is a en- 
ergy scale separation between the quasi-particle and the 
vortices: the vortices are moving much slower than the 
quasi-particles, therefore can be treated in a hydrody- 
namic approximation. The important effects of Volovik 
and AB phase scattering on ARPES data was discussed 
in Ref. jSj. 

As explicitly stressed in this paper, the U(l) mutual 
CS theory in the continuum does not have the period- 
icity under 9 — > 9 + 1, because moving hc/2e vortex is 
different from moving —hc/2e vortex. In fact, the peri- 
odicity of the conventional CS theory under 9 — ► 6 + 2 is 
also a very intricate issue. On the one hand, the pertur- 
bative RG expansion in terms of the statistical angle 9 
in the continuum limit in Refs. |^ , j3l|]32[ does not have 
this periodicity. On the other hand, properly regularized 
on the lattice, the CS theory does have this periodicity 
p3| . The two different regularizations do lead to differ- 
ent conclusions on the Quantum Hall transitions. It also 
leads to notorious Hall conductivity difficulty at v = 1/2 
Ej . A deep understanding on the periodicity problem in 
both conventional CS term and mutual CS term is clearly 
required in the near future. 
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APPENDIX A: THE MOST GENERAL FORM 
OF THE GAUGE PROPAGATORS 

Adding the gauge fixing terms ^{(d^af) 2 + (<9 M a*) 2 ) 
to Eqn.^5|, we can find the gauge field propagators by 
inverting the matrix 



ni(fc)M<W — W~ > 



n 3 (fc)fc(<w -b^L ) + b^ 

(Al) 



The results are: 
!r< \ — a -k il K v 

[S'^i^u — —^4 r 

(G<jxj>)^i/ 



n. 



1 / c k^k v , 



n 2 k 



k 2 



n, 



n 1 n 3 + n2 ( 



ni 



k 4 



nin 3 



k 2 ' 

1 / c k^k v . 



n 2 k 



k 2 



In the Landau gauge a = and putting k 
recover the results calculated in Sec. III. 



(A2) 
0, we 



APPENDIX B: MOVING HC/E VORTICES AND 
THE IRRELEVANCE OF THE DYNAMIC 
DOPPLER SHIFT EFFECT 

For hc/e vortex, we can perform a neutral-like single- 
valued gauge transformation by setting </>a — = 0/2. 
The double strength vortex and its stability was inves- 
tigated in Ref. the associated spin-charge separation 
has been discussed extensively in Refs. |l],[Hj. In this 
appendix, considering a moving hc/e may still scatter 
quasi-particles by the AB phase, in analogy to the hc/2e 
vortices discussed in the main text, we still include a pos- 
sible mutual C-S interaction with the statistical angle 9 
taking any integer values. For 9 = 0, we recover the 
results in Ref. 0] and also get the new result on the ex- 
pression of current operator 1^ in the superconducting 
phase. Although the static Volovik effect is marginally 
relevant [|6po[| and cause finite density of states of the 
quasi-particles at zero energy. In this appendix, we show 
that surprisingly the dynamic Volovik effect is highly ir- 
relevant on the both sides of the QCP. 

Putting 4>a = 4>b = 0/2, namely setting a M = in 
Eqn. 
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leads to: 



£ = V'lal^T + VfP x T 3 + VAPyT^lpU + (1 -> 2, X 



K 



(Bl) 
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where the effective gauge field is Afff — — K~ l Jfj,. 
There is no internal gauge field a M , because static double 
strength vortices do not scatter the quasi-particles by AB 
phase. 

Performing duality transformation to vortex represen- 
tation leads to 



C = i>hM a + - ia^)<P\ 2 + F(|$|) 



(B2) 



where the last ( Berry phase ) term can be absorbed into 
A e Jf by redefining Afff — > Afff + i/LtJ M o- Comparing 
with Eqn.^2], there is only charge fluctuation. There is 
no mutual C-S term, namely 9 = 0. 

The duality transformation to boson representation is 

c = ^hM a + \(d, - ia^\ 2 + v(\4>\) 

+ \(fp 2 + ia^M4 ~ A x f ) + \fl (B3) 
Integrating out the electric gauge field leads to 



(B4) 



Setting = Aff* up to irrelevant term leads to 

C = Tph^a + fltfdrt +\{d M - lA^\ 2 + V(\cb\) + ■ ■ ■ 

(B5) 

where the second term is the Berry phase term which 
is a first order time derivative and • • • means the irrele- 



vant couplings between bosons and fcrmions. Eqn.B5 is 



essentially the same as Eqn. Bl 



However, a moving hc/e vortex may still scatter quasi- 
particles by AB phase. In general, we should also in- 
clude a possible mutual C-S term in the effective action. 
The most general form is again given by Eqn.|2^ with 
the statistical angle 9 taking any integers and Aff' = 
A^ — K~ 1 J fi for hc/e vortex. The formulation developed 
in Sec. IV can be straight- forwardly applied to hc/e vor- 
tices. Putting the above two values for 9 and A e Jt into 
Eqn|5^ leads to 



- iOA^ipa 



(B6) 



where the second and third terms are the Berry phase 
terms which are first order time derivatives. 

The above equation indicates that the quasi-particlc 
carries charge 9 which can take any integer numbers. 
Obviously, there is no periodicity in 9. The most nat- 
ural and simplest choice is 9 = where the above equa- 
tion reduces to Eqn.BS. In this equation, bosons and 
fermions are asymptotically decoupled. The fermions 



(spinons) carry only spin 1/2, the bosons ( holons) carry 
only charge e. There is spin-charge separation. The 
long-range Coulomb interaction could be incorporated by 
adding the dynamic term for the time component of the 
gauge field ^fc|Ao(fc)| 2 Ej|. Without the fermionic part, 
the action is the same as the superconductor to insula- 
tor transitions studied in p6[ , except here condensed is 
a charge e boson instead of a Cooper pair. 

In the underdoped regime, the vortex condenses < 
<f> >= $0 which generates a mass term for a M in Eqn.B2 
i^o( a A»)f- Integrating out the massive a M leads to a 



2 ^oV"m)* • 

Maxwell term for A e J f = A„ - K~ 
can identify the automatically conserved electric current 
/„ = ■Jvid^Jfi - d^dvJv) ||. The dynamic Dopplcr 



1 Jf Jl from which we 



shift effect is ~ ~i^(dJ) 2 , it is a four fermion term with 



two derivative, therefore highly irrelevant |jj . The static 
Volovik effect is marginally relevant ||[l0| and cause fi- 
nite density of states at zero energy. Surprisingly, the 
dynamic Volovik effect is highly irrelevant. Note that its 
coefficient $^~ 2 diverges as we approach the QCP from 
the underdoped side. 

In the superconductor phase, the vortex <!> is mas- 
sive, therefore can be integrated out, it leads to the 
old Maxwell term j^(f^) 2 - Integrating out gen- 
erates a mass term ^m^{A e J^) 2 . From the mass term, 
we can identify the automatically conserved electric cur- 
rent I^i = m$( Jfi)t- Note that although the quasi-particle 
electric current is not conserved itself, the total electric 
current 1^ = m$((5 Jl(I / — ■^r L )Jv{q) is conserved. The 
dynamic Doppler shift effect due to the virtually fluctu- 
ating vortex-anti- vortex pairs is ~ to$(J m ) 2 . It is a four 
fermion interaction term, therefore irrelevant. Note that 
its coefficient vanishes as we approach the QCP from the 
superconducting side. 

We reach the conclusion that the dynamic Doppler 
shift term is irrelevant on both sides of the QCP. We 
conjecture that it is likely also irrelevant at the QCP. 
Counter-intuitively, it contains two more derivatives in 
the vortex condensed underdoped regime than in the vor- 
tex depleted superconducting phase, therefore even more 
irrelevant in the underdoped regime than in the super- 
conducting phase. The previous scaling argument on the 
irrelevance of this term at the QCP presented in Ref. 
]jj is questionable, because the authors scaled the phase 
(f) like a free field. In fact, <j> should be confined to be 
< (j> < 2ir, so it is far from being a free field. 

Putting (j)A — 4>b = 4>/2 m Eqnj2^, we get the electron 
annihilation operator: 



e a pe^ 2 e~^ 



(B7) 



The above Eqn. demonstrates clearly that electron is 
separated into charge e boson b = e l ^l 2 and spin 1/2 
spinon ip. Note that the crucial difference from Eqn.|25| 
is that there is no internal gauge field ! 
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The electron Green function is given by: 

G(x,t) =< b(x,tp (0,0) > [e lR >- 3 < Vii a (^*)V>L(0,0) 



-iKi -x 



<^\ 2a (S,t)^ i2a (0,0) >] 



In the underdoped regime, boson b = d"**' is massive, 
therefore short-ranged correlated. Being the product of 
fermion and boson Green functions, the electron corre- 
lation function is also short-ranged. The underdoped 
regime is a Mott insulator as suggested in Refs. ]i], |i"9| , p0[ . 

At the critical point, in real space, b decays with a 
power law 1 + rjxy, the electron Green function decays 
with a power 2 + 1 + tjxy ■ 

In the superconductivity phase, < b >=/= 0, the electron 
Green function decays with a power 2. The anomalous 
Green functions F(x, t) =< C-\ (x, t)Ci (0, 0) > also starts 
to form: 
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